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Abstract 


When the predecoupling plasma is thermodynamically reversible its fluctuations are clas¬ 
sified in terms of the adiabatic and entropic modes. A different category of physical solutions, 
so far unexplored, arises when the inhomogeneities of the viscosity coefficients induce com¬ 
putable curvature perturbations. The viscous modes are explicitly illustrated and compared 
with the conventional isocurvature solutions. 
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The adiabatic and the entropic modes have been introduced long ago [Tj (see also [2] 
and references therein) in connection with the pioneering analyses of the temperature and 
polarization anisotropies of the Cosmic Microwave Background (CMB in what follows). The 
relative position of the first anticorrelation peak of the cross-spectrum between the tem¬ 
perature and the polarization demonstrated (already from the first releases of the WMAP 
data that the initial conditions of the Einstein-Boltzmann hierarchy are predominantly 
adiabatic. This evidence has been further scrutinized and quantitatively refined by the 
subsequent WMAP releases and by Planck explorer results |4]. While the conventional dis¬ 
tinction between the adiabatic and the entropic solutions posits that the ambient fluid is 
thermodynamically reversible, it seems plausible to relax this hypothesis with the aim of 
complementing and extending the physical initial conditions of the predecoupling plasma. 
According to the current data M, anticorrelated non-adiabatic modes in the presence of 
a dominant adiabatic mode may even improve the fit of the temperature autocorrelations 
accounting for potential large-scale suppressions of the corresponding angular power spectra. 
This same strategy could also be investigated in the case of the viscous mode derived in this 
analysis. 

The fluctuations of the total pressure of the predecoupling plasma are customarily decom¬ 
posed into an adiabatic component supplemented by the entropic (or simply non-adiabatic) 
contributions: 


5p t = c 2 st 5p t + 6p nad , 


""St 




( 1 ) 


where the subscripts in the brackets signify that the variations must be performed, respec¬ 
tively, at constant specific entropy and constant energy density pt . If the cold dark matter 
(CDM in what follows) is only supplemented by a radiation background, the specific en¬ 
tropy is the ratio between the entropy density of the radiation and the CDM concentration, 
i.e. approximately ~ T 3 /n c . The entropy contrast S = can then be written as 
S = (36 r /4 — 5 C ) where 5 r and 5 C are, respectively, the density contrasts of the radiation and 
of the CDM species. If S = 0 the density contrasts are in a fixed ratio (i.e. S r /5 C = 4/3) and 
the solution is adiabatic; if S ^ 0 the adiabaticity condition is violated and 5 r ^ 4<5 c /3 (see 
e-g- 0). Across the matter-radiation transition the entropic solutions are determined by the 
total sound speed of the plasms^ and by the specific form of 6p nad ; for the CDM-radiation 
mode we have: 

A 

SPnad = PcC 2 st S, ( 2 ) 


c 2 = 

c st 


3(3a + 4) ’ 


where p c denotes the CDM energy density a = a/a eq is the scale factor normalized at the 
matter-radiation equality. The analytic expression of 5p nad can be generalized to a generic 
pair of species of the plasma (see e.g. fourth article in Ref. 0)- 


2 We shall consider hereafter the case of a conformally flat background geometry characterized by a scale 

factor a(r ) where r denotes the conformal time coordinate. 
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Even if the conventional terminology might suggest otherwise, the non-adiabatic modes 
arise in a globally inviscid fluid, as the preceding considerations illustrate. In the present 
paper we want to drop this hypothesis since the total energy-momentum tensor of the plasma 
could include the contributions of the shear viscosity, of the bulk viscosity and of the heat 
transfer. The adiabatic limit, in a strict sense, is recovered when the viscous contributions 
are neglected and the total entropy four-vector is conserved. Across the matter-radiation 
transition the shear viscosity coefficient rj determines the optical depth, the Silk damping 
scale and, ultimately, the shape of the visibility function met Conversely the role of the 
bulk viscosity coefficient £ is disregarded since it must vanish in the limit c s t —> l/\/3 [6j. 
The scaling properties of the viscous coefficient can be expressed, in a unified notation, as: 


£ £ /I 

_ I _ 

V~ T X~ \ 3 




(3) 


where x denotes the coefficient of heat conduction; the values of q > 1 parametrize the rate 
of extinction of the bulk viscosity coefficient prior to matter-radiation equality, i.e. in the 
limit r r eq . It is possible to argue [6], with reasonable arguments, that q ~ 1 in Eq. ((3]) ; 
however, to avoid supplementary hypotheses, the values q will be taken as a free parameter. 

Prior to matter radiation equality but well below the temperature of neutrino decoupling, 
the total energy-momentum tensor of the fluid shall then be written, for the present ends, 
as: 

V (p t , n t , £) = (p t + p t )u^u v - p t 8" + £ (4) 

The total pressure of the plasma pt is in general defined to be the same function of p t (i.e. 
the total energy density) and nt (he. the total concentration) as in the adiabatic limit. 
The Eckart approach^ seems preferable in the present situation where the concentration 
of radiation quanta exceeds the concentrations of the other species. The first principle of 
thermodynamics together with the covariant conservation of the total energy momentum 
tensor and of the global particle current (i.e. j a = nt u a ) implies the law of growth of the 
entropy and the rate of suppression of the temperature: 


= |(V Q u“) 2 , 


T = -Tc s 2 t VX, 


(5) 


where S' 1 = n t S u M is the entropy four-vector in the case when the bulk viscosity is the only 
source of irreversibility^. 

The bulk viscosity coefficient £ is the sum of a homogeneous part supplemented by the 
inhomogeneous contribution: 


£(f,r) = a(r)£(r) + a(r)<5£(f, r). (6) 

3 In the Landau approach the four-velocity coincides with the velocity of the energy transport while in 
the framework of the Eckart theory denotes the velocity of the particle transport [7]. 

4 Note that T is the temperature of the plasma that coincides, roughly speaking, with the photon tem¬ 
perature since the concentration of the photons greatly exceeds the concentration of the other species both 
relativistic and non-relativistic. 
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The evolution equations for the homogeneous expansion rate are a consequence of Eq. (HI) 
and of the Einstein equations: 

(V 2 - %') = 4vrGa 2 (pt + Pt), 3 U 2 = 8vr Ga 2 p t , (7) 


where Vt = Pt — 3£ "H is the shifted background pressure. The gauge-invariant bulk viscosity 
fluctuations are given b}@: 

Z = 5Z+[t + m\(B-E'), ( 8 ) 

where B and E parametrize the off-diagonal (scalar) degrees of freedom of the perturbed 
metric. If £(r) = £ (t) = 0 then 5 £ equals 5 and the bulk viscous fluctuations are automati¬ 
cally gauge-invariant; in the same interesting limit Vt — > Pt and the shifted pressure coincides 
with the pressure itself. The Hamiltonian and the momentum constraints for the evolution of 
the fluctuations are not modified by 5. On the contrary 5 appears in the ( ij) components of 
the perturbed Einstein equations. Consequently the evolution of the curvature perturbations 
on comoving orthogonal hypersurfaces is given by: 


<g + m K + 9) 


m 2 „ ^ h 

Pt + V t ^ + P + V t 


Ucl 

4nGa 2 (p t + V t ) 


V 2 T, 


(9) 


where <f> and T denote the Bardeen potentials entering the definition of the curvature pertur¬ 
bations as 77 = -[T + P(P4> + T')/(P 2 ■ InEq. ® 0 t = V • v t is the three-divergence 

of the gauge-invariant velocity field. On the right-hand side of Eq. (j9]) the only term surviv¬ 
ing the large-scale limit is the first one; V 2v k is sub leading and, similarly, © t can be neglected 
thanks to the momentum constraint implying that 

yvy+w$) _ w(k+£) 

01 (U 2 -W) H ( 0) 


where the second equality in Eq. (fTO follows from the gauge-invariant definition of 1Z. 
Equation (1TU|) shows, as it is well known, that the total velocity field can be neglected in the 
large-scale limit. 

The total viscosity is given by S(fc, r) = S 7 (/c, r) + H„(A;, r) and it is proportional to the 
total radiation energy density pn = p 1 +p u - 1 note that S 7 (fc, r) and S v (fc, r) are the viscosities 
of the photon and of the neutrino sectors. According to Eq. fl3J) S 7 and go to zero deep 
in the radiation epoch and they must be proportional to the corresponding energy densities; 
thus they can be parametrized, in Fourier space, as: 


h(k,T) = ^ 


n [ 


3 - 4 (~) 


2g 


V*(fc), 


'(k, T ) = 


3 - c st( r ) 


2 q 


v*(fc), 


( 11 ) 


5 We decompose the scalar fluctuations of the four-dimensional metric tensor g as Sgoi = —a 2 diB 1 
5gij = 2a 2 {tp5ij — didjE) and Sgoo = 2a 2 <j>. The two gauge-invariant Bardeen combinations correspond to 
$> = (/)+ (B — E')' + U{B - E') and T = tp - H(B - E'). 
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where V*(/c) denotes the fully inhomogeneous contributiorj§. 

Barring for a potentially important role of the bulk viscosity coefficient in the very early 
Universe, as repeatedly suggested in the past (see for instance M), £ cannot appreciably 
affect the evolution of the homogeneous expansion rate across the matter-radiation transition. 
We shall therefore accept that £ = = 0, so that the solution of Eq. (J7D across the matter- 
radiation transition will be given, in conformal time, by a = a/a eq = ( x 2 + 2x) where 
x = t/t\ and r eq = (y/2 — 1 )t\. Equation (JUJ) can then be solved in the large-scale limit 
using a as integration variable and recalling that, by definition, T~L = din a/da] the result 
of this step can be written as: 


9 


a 


2 q 


K(t,a)=7J.(t) + 5 jV.(fc). 


( 12 ) 


In Eq. (fT2|) 77*(fc) denotes the conventional adiabatic mode while V*(k) corresponds to the 
viscous mode. Recalling the definition of 77 in terms of T, the expression for d/( k , a) can be 
derived in terms of hypergeometric functions and the two relevant asymptotic limits (well 
before and well after equality) are: 


lim T(/c, a) 

a<Cl 


lim ddfc, a) 

a> i 


->■ 




-~Kjk) 

~n*(k) 


1 ■ d +o( “ 2) 
1+ Ju°(“") 
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4q(2q + 3) 

q 3-2 q 


V.(k) 


Ot\ 

4/ 


2g 


( 13 ) 


The derivation of Eq. (1T3]) assumes that = T and implies the absence of any source 
of anisotropic stress; while this approximation holds over large-length scales, the presence 
of the neutrinos introduces a difference between and T. The full structure of the viscous 
mode will now be obtained in the limit a 1 and in the presence of the neutrinos, taken to 
be massless as in the context of the concordance paradigm. The gauge-invariant evolution 
equations for the lowest multipoles of the neutrino hierarchy in the presence of viscous 
inhomogeneities can be written as: 


S' = W + 


9 H 2 
Pv 


- -0,. 


0' v = ^—8 V - k 2 cr u + k 2 & - k 2 : 

4 4 p v 


cr. 


= — 0 V - 
15 


10 ' 


-Tr 


u3i 


( 14 ) 

( 15 ) 


where T v 3 is the octupole of the perturbed phase space distribution of the neutrinos. In the 
tight-coupling limit the evolution equations of the baryon-photon system are instead: 

/2 4 


gv 2 

s r = w + —s 7 - 

7 p 7 

UR b rj 

-V~fb + 


- 9 . 


7 bi 


Q'-yb + 


k 2 9 . 


7 b 


Rb + 1 


Py ( Rb + 1 ) 


<51 = 3T' 


= k 2 $ + 


9'H 2 

6 lb + —± b , 


Pb 


k 2 S „ 


9 n 


4 (Rb + 1 ) 4 p b 


k 2 (^b + “ 7 )- 


( 16 ) 

( 17 ) 


6 If the inhomogeneities are localized either in CDM or in the baryon sector Eq. CHt) can be complemented, 
respectively, by S c and hereafter we shall focus, for illustrative purposes, on the case where Sf, = S c = 0. 
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In Eq. (TTTp . as usual, R b = 3pb/(4p 7 ) denotes the ratio between the energy density of the 
baryons and the one of the photons; the shear viscosity coefficient. Finally the equations for 
the CDM fluctuations are simply 

o' c + ne c = k 2 $, s' c = 3d ,/ - e c . (is) 


We are looking for a solution of the system formed by Eqs. (J9]) , (jT4|l — (1T5|) . (TT 6 |l — (Tl7|) and 
(TT 8 |) in the region kr kri <C 1 (and r < Ti) that corresponds to wavelengths larger than 
the Hubble radius prior to equality. The Bardeen potentials become, in this double limit 
^(kpx) = T if {k)x 2q and <h(fc,a;) = <& Jf (k)x 2q and the corresponding amplitudes are related 
V*(fc) as: 

q 2 2q+2 [$*(k) + 2 (q + 1 )**{k)} + 9V*(A;) = 0. (19) 

Furthermore, the direct integration of Eqs. (jT4T) . (fl 6 |) and (1T8D leads to the following results: 


5 7 (fc, x) = 5 u (k, x) ~ 


4 T*(fc) + 


q 2 2q+1 


V,(fc) 


x 2q , S c (k,x ) = S b (k,x) — 3 ^*(k)x 2q . ( 20 ) 


Inserting Eq. (120D into Eqs. (fT5j) and (1171) the three-divergences of the peculiar velocities 
can be obtained: 


0 u {k,x) = 


k 2 n 


( 2 g+l) 


0 lb (k,x ) = 


k 2 Ti 


4/*(fc) 4- $ Jt (fc) + 


4t*(^) T *&*(&) + 


(2<Z+1) 

t .2 

0c(k,x) = 2 ( q+ i) x2q+1 - 


9(1 ~ 2 q) 

q 22g+3 

9 

q 2 2 | 3+ 3 


V*(k) 

3 q 


x 


2q+l 


1 - 


2 R b 


V*(fc) 


x 


2q+l 


( 21 ) 


Finally direct integration of Eq. (ITS]) implies an explicit expression for the dimensionless 
anisotropic stress of the neutrinos: 


cr u (k,x) 


2fc 2 7f[ff»(fc) + $«( fc )] 2( q +i) 3 k 2 r 2 n - 2 q 

15(g +l)(2g +1) X 5 q 2 2q + 2 V1 + 2q 


V,(fc) x 2(q+1 \ 


( 22 ) 


where the hierarchy of the neutrinos has been truncated by requiring T v 3 = 0, as it happens 
in the case of the adiabatic mode and of some of the entropic solutions. Equation (|22l) 
must also be compatible with the (i =4 j) component of the perturbed Einstein equation, 
i.e. V 4 ($ — T) = 12nGa 2 d$jWi\ in the concordance paradigm the total anisotropic stress 
coincides with the neutrino component and it is given by didjH 1 -* = Ap u X7 2 cr u /3. This relation 
generates a further condition on a v i.e. 6R u a u (k,x) = /c 2 rf(T* — <h*)a; 2( '' J+1 ' > where R u = 
p v / (p 7 + p u ) denotes the neutrino fraction in the radiation plasma. Using this expression for 
a u together with Eq. (122j) it is easy to obtain an algebraic relation between <h*(fc), 4/*(/c) 
and V*(/c); such a relation can be inserted back into Eq. (HUD to derive <h*(fc) and 4/* (k) in 
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terms of V*(/c): 


9V,(fc) [q {AR U — 5) — 5] 

2 %+!)^ [4 J R I/ + 5(g + l)(2g + 3)]’ 

9 V*(fc) [2i2„ + 5(g + 1 )] 

2 %+Dg [4R i , + 5(2g + 3)(g + l)]‘ 


(23) 


For the sake of accuracy it is finally appropriate to mention that the solution derived here 
satisfies the Hamiltonian and the momentum constraints stemming fro the (00) and (Of) 
components of the perturbed Einstein equations. 

The analytic form of the viscous mode for typical wavelengths larger than the Hubble 
radius prior to matter-radiation equality can be used as initial condition of the Einstein- 
Boltzmann hierarchy. To illustrate this point and to spell out more clearly this strategy the 
power spectra of the adiabatic mode, of the viscous mode and of their cross-correlation can be 
assigned at the common pivot scale k p = 0.002 Mpc -1 which corresponds to £ p ~ 30 in multi¬ 
pole space. The power spectra will be denoted by the unified notation Px(k) = Ax{k / k p ) ny ~ l 
where X = TZ, V, 1ZV are, respectively the amplitudes of the adiabatic mode, of the viscous 
mode and of their mutual correlation; similarly y — r, v, rv denote the corresponding spectral 
indices. The temperature anisotropies for multipoles £ < ^/z rec (where z iec ~ 1090.51 is the 
redshift of recombination) are unaffected by the thickness of the visibility function |5j. The 
SW contribution typically peaks for comoving wavenumbers k ~ 0.0002 Mpc -1 while the inte¬ 
grated Sachs-Wolfe effects contributes between k min = 0.001 Mpc -1 and fc max = 0.01 Mpc -1 . 
Recalling that 5 1 = 4A/o where A/o is the monopole of the intensity of the radiation held, 
from Eq. (fTTfl) we have that 


Sy(k, Tree) 5ry(k, Tj) 




(24) 


where r* corresponds to the initial time at which the normalization is set (typically Tj -C r eq 
and cq <C 1). The SW contribution is then given by: 


A[ SW) 


(fc,/i,T 0 ) = 

/(<?) = 


+ f{q)v*{k) 


= — ifikr 0 


S 3-2 q 

20 q 


^ 5 / 3ct rec 

_ 6 V 4 


2 q 


F[2q, 2 q, 1 + 2 q, -3a rec /4] L 


(25) 


where F[a, b, c, z] denotes the hypergeometric function and fi = k ■ h is the projection of 
the direction of the Fourier mode in the direction of the photon momentumQ Following 
a standard procedure and denoting k 0 = r$ l the angular power spectrum of the full SW 
contribution will be given by: 

CT } = £) + f 2 (q)A v Z v (n v , £) + AnvZnv{n rv , £) cos 7 ^, (26) 

25 5 

7 In the limit a re c —> 0 the second term in the bracket goes to zero. In the limit a rec 1 the same term 
provides a logarithmic correction. To simplify Eq. (12511 the limit cu = 0 has been assumed. 
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T- 


(27) 


Zx{n x ,£) 


7T 2 /k 0 \ nx 1 

T Vv 


r(3- n jr(<+f -i) 

r 2 ( 4 ¥*)r(l + <-^)’ 


where 7 r „ denotes the correlation angle. Equations (|26l) and (127)1 can be used as starting point 
for dedicated analyses where the dominant adiabatic mode is correlated (or anticorrelated) 
with the viscous mode. 

The viscous curvature perturbations vanish faster than in the case of the conventional 
entropic modes. The baryon-radiation and the CDM-radiation modes lead to Bardeen po¬ 
tentials and curvature perturbations that vanish as a in the limit a « 1; in the same limit 
the curvature perturbations induced by the viscous mode vanish as o? q with q > 1. The 
curvature perturbations corresponding to the neutrino entropy mode are identically zero 
(to leading order) while the Bardeen potentials are approximately constant. The neutrino 
velocity mode arises when the momentum constraint is vanishing at early times and, simul¬ 
taneously, the total density is uniform. This requirement, once inserted into the evolution 
equations in the longitudinal gauge, entails a singularity in the Bardeen potentials; this 
mode (singular in the longitudinal description) is regular in the synchronous gauge and has 
no direct counterpart in the viscous case. These considerations show that prior to matter- 
radiation equality the viscous solution is really and truly of isocurvature type since, in this 
limit, the curvature perturbations tend to zero faster than in the case of their (regular) 
entropic counterparts. 

All in all the results of this investigation demonstrate the existence of a regular class of 
viscous solutions for the fluctuations modes of the predecoupling plasma. Since the CMB 
anisotropies are sensitive not only to the homogeneous matter content of the Universe but 
also to the Cauchy data of the cosmological perturbations, the viscous modes could serve as 
sound initial conditions of the Einstein Boltzmann hierarchy in conjunction with a dominant 
adiabatic component. In the present framework the direct observational limits stemming 
from the temperature and polarization anisotropies of the CMB can be translated into specific 
constraints on the properties of the bulk viscous stresses across matter-radiation equality. 
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